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Abstract 

MINOS is a long baseline neutrino oscillation experiment that uses two de- 
tectors separated by 734 km. The readout systems used for the two detectors 
are different and have to be independently calibrated. To verify and make a 
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direct comparison of the calibrated response of the two readout systems, test 
beam data were acquired using a smaller calibration detector. This detector 
was simultaneously instrumented with both readout systems and exposed to 
the CERN PS T7 test beam. Differences in the calibrated response of the 
two systems are shown to arise from differences in response non-linearity, 
photomultiplier tube crosstalk, and threshold effects at the few percent level. 
These differences are reproduced by the Monte Carlo (MC) simulation to 
better than 1% and a scheme that corrects for these differences by calibrat- 
ing the MC to match the data in each detector separately is presented. The 
overall difference in calorimetric response between the two readout systems 
is shown to be consistent with zero to a precision of 1.3% in data and 0.3% 
in MC with no significant energy dependence. 

Key words: neutrino detector calibration, iron-scintillator sampling 
calorimeter, test beam measurements, readout system 
PACS: 29.40. Vj, 29.40.Mc, 29.40. Gx 



1. Introduction 

The Main Injector Neutrino Oscillation Search (MINOS) is a long base- 
line, two-detector neutrino oscillation experiment that uses the NuMI neu- 
trino beam at Fermilab jH . The energy spectrum and flavour composition of 
the neutrino beam is measured at two detectors located 734 km apart: the 
Near detector (ND) at Fermilab and the Far detector (FD) at the Soudan 
Underground Laboratory in Minnesota j^. MINOS is sensitive to neutrino 
oscillations in the region studied by the atmospheric neutrino experiments 
and has recently measured | Am^| = 2.43 ± 0.13 x 10"^ eV^ jif. 

To control systematic errors, such as those that arise from uncertainties 
in neutrino flux and cross-sections, MINOS was designed as a two-detector 
experiment. This design allows a relative measurement to be made but it 
subsequently becomes necessary to ensure a precise inter-detector calibration 
of the energy scale. The MINOS ND and FD have a similar steel-scintillator 
structure but were instrumented with different readout systems. The driving 
factors for using different readout systems were the relatively small size and 
high event rate of the ND compared to the large size and low event rate of 
the FD. To investigate the effects of the two readout systems on the relative 
energy scale, test beam data were taken with the dedicated MINOS Calibra- 
tion detector (CalDet) [4] that was simultaneously instrumented with both 
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readout systems. Specifically, differences arising from response non-linearity, 
photomultiplier tube crosstalk, and threshold effects were studied. 

This paper quantifies the differences between the two readout systems 
and demonstrates how they can be controlled. In Sections [2] and [3] the ex- 
perimental setup and the selection of the events used in the analysis are 
described. Section H] outlines the calibration scheme and Section O details 
the uncertainties on the calibration. Comparisons of strip occupancy and 
PMT crosstalk for the two readout systems are given in Section [61 The re- 
sults of the calorimetric response comparisons and the conclusions are given 
in Section [7] and Section [8] respectively. Further details of this analysis can 
be found in [sj. 

2. The MINOS Calibration Detector 

The MINOS detectors are tracking-sampling calorimeters consisting of 
alternating layers of steel and scintillator planes. CalDet planes use 2.50 cm 
thick unmagnetised steel in contrast with the ND and FD that use 2.54 cm 
thick magnetised steel. In all three MINOS detectors, the plastic scintilla- 
tor planes are made out of strips which are 1 cm thick and 4.1 cm wide. 
Successive scintillator planes are rotated by 90° to allow three dimensional 
event reconstruction. CalDet consists of 60 steel planes measuring 1 x 1 m 
square, interleaved with scintillator planes comprising of 24 strips of 1 m 
length. Light produced in the scintillator is captured by wavelength shifting 
(WLS) fibre optic cables. At two corners of each CalDet plane the WLS 
fibres from the 24 strips are brought together in a manifold. Additional fibre 
optic cables are then attached at the end of the manifolds to guide the light 
to multi-anode photomultiplier tubes (PMTs). The light output of CalDet 
was engineered to match the ND and FD by using long fibres: 6 m clear 
fibres on one end of the strips and 3 m green WLS fibres on the opposite 
end. 

A schematic diagram showing the configuration of CalDet when it was 
simultaneously instrumented with both readout systems is shown in Fig. [TJ 
The view of the detector shown is looking downstream at the first plane with 
successive planes following into the page. Each scintillator strip was read out 
at both ends: one end was instrumented with the ND readout system while 
the opposite end used the FD readout system. Successive planes alternate 
between the ND readout system being connected via clear fibres and the FD 
readout system by green WLS fibres and vice versa. 
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Figure 1: A schematic diagram showing the first 6 planes of CalDet. In the configuration 
shown, CalDet was simultaneously instrumented with both ND and FD readout systems. 
Clear fibres are connected to the two upper manifolds and green WLS fibres to the two 
lower manifolds. 

The ND uses 64-anode Hamamatsu M64 PMTs Q with a custom base 
divider ratio of 3-2-2-1-1-1-1-1-1-1-1-2-5 and electronics based on the QIE 
ASIC [?]. The FD readout system uses 16-anode Hamamatsu M16 PMTs jsf 
with a custom base divider ratio of 2.4-2.4-2.4-1-1-1-1-1-1-1-1-1.2-2.4 and 
electronics based on the VA ASIC ^]. 

The multi-anode PMTs used in MINOS have a single photocathode. The 
region of the cathode from which charge is focused on to a specific anode is 
called a pixel. The PMTs were housed in custom made boxes with optical 
routing that guided the light from individual strips to specific pixels. The 
routing pattern was designed such that neighbouring strips did not use neigh- 
bouring pixels. This was done to minimise the effect of crosstalk on pattern 
recognition. Crosstalk predominantly occurs on pixels next to the one illu- 
minated due to leakage of photoelectrons before the first dynode [gI, Isj. The 
optical routing of the FD readout system was such that three PMTs were 
fully utilised to readout two planes. The ND PMTs have 64 pixels, hence 
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only one PMT was required to read out two detector planes and 16 pixels 
remained unilluminated. 

In the ND front-end electronics, the charge from each PMT pixel is inte- 
grated and digitised continuously at 53 MHz (every 18.8 ns) for up to several 
tens of /IS with no dead-time. The digitisation is performed by a specialised 
ASIC, the QIE-chip. It splits the incoming signal into eight binary-weighted 
"ranges" and integrates the resulting fractional currents on capacitors in each 
range. A sample-and-hold circuit stores the integrated charge in each range. 
The QIE-chip then selects the first non-saturated range and passes its held 
voltage to the 8-bit ADC for digitisation. This combination of the 8 ranges 
with the 8-bit ADC gives the QIE-chip the same dynamic range as a 16-bit 
ADC but with an approximately constant fractional error across the entire 
dynamic range. 

The QIE circuit does not provide a linear response by itself because each 
range has a different ADC offset and gain. A 16-bit DC current injection 
circuit is used to map the response of the chip. The results of that cahbration 
are uploaded locally in a look-up-table and used during normal data taking. 
The calibrated output of the ND electronics is linear over the entire dynamic 
range to better than 0.5%. The stored data are processed once the acquisition 
period is complete. Pedestal subtraction and linearisation are performed for 
each channel via the look-up-table in a single step. 

The FD readout system is based on the front-end ASIC VA32_HDR11 
(short VA chip), developed in collaboration with the Norwegian company 
IDEA ASA. The chip includes a charge sensitive preamplifier, a shaper, sam- 
ple and hold for each channel, and is followed by an analogue output multi- 
plexer. The output from a selected channel can be switched to a differential 
output buffer, which drives a significant length of cable to a remote ADC. 
The ADC value is proportional to the integrated charge of the PMT pulse, 
but starts to saturate between 15-30 pC. Readout of the front-end electron- 
ics is triggered by a signal from the last dynode of each PMT and the dead 
time after each digitisation is 6 //s. This electronics was a cost effective choice 
for the low rate environment of a deep underground detector. 

The amount of recorded data is reduced in both systems by recording 
hits above a readout threshold of about 1/3 of a photoelectron (PE). For the 
FD readout, this threshold is apphed to the total integrated charge, while 
for the ND it is apphed separately to each 18.8 ns integration sample. 
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2.1. The PS Test Beam and External Trigger 

The results presented in this paper are based on data obtained by ex- 
posing CalDet in the T7 test beam in the East Experimental Hall of the 
24 GeV/c CERN Proton Synchrotron (PS) The mixed composition 

beam of e, fj,, n and p was operated between 0.6- 10 GeV/c. The beam line 
was instrumented with several Cherenkov counters to identify electrons and 
a time-of-flight (TOP) system to further aid particle identification. Control 
of the beam momenta spread (Ap/p < 2%) and the instantaneous event rate 
(< 1 kHz) was achieved by adjusting the beam line collimators. Por further 



details see m, 11 



Synchronous readout of both systems was accomplished by an external 
trigger provided by the TOP counters. After a trigger, the ND and PD 
electronics were enabled for 376 ns and 1 /is time windows respectively, which 
was long enough to record individual events. In addition, for 50/is after each 
trigger, signals from the TOP system were suppressed to ensure that the 
electronics had ample time to complete the digitisation of the first event. 

2.2. The Monte Carlo Simulation 

The simulation results presented alongside the data in this paper use the 



MINOS Monte Carlo (MC) simulation. The MC is based on GEANT3 [12 
and is used to generate raw energy depositions [GEANT hits), which serve 
as the input to a C++ based detector response model. The response model 
includes the effects of light collection and propagation from the scintillator 
strips to the PMT photocathode as well as the PMT and the electronics. 

The attenuation of light along the length of a strip is modelled as the 
sum of two exponential functions to give a short and a long attenuation 
component. Pits to both bench-top and cosmic ray data determined the 
lengths and relative contributions of the long and short components of the 
attenuation. In the simulation of the CalDet the approximation is made 
that every strip has the same short and long attenuation components. The 
variation in the light output of the strips and the gains of the PMTs are 
directly incorporated into the simulation by using the calibration constants 
determined by the detector calibration procedure (see Section H]). 

3. Event Selection 

The readout system comparison required separation of specific particle 
types and single-particle interactions in CalDet. Cuts to achieve this were 
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developed and are presented here. Positron and muon samples were chosen 
because they are efficiently identified with high purity and their energy loss 
mechanisms are well understood. When passing through iron at GeV-scale 
energies the positrons form a dense electromagnetic shower, which makes 
them ideal candidates for probing the full dynamic range of energy depo- 
sitions expected for neutrino interactions in MINOS. The muons were used 
to probe some of the systematic errors of the calibration, their energy loss 
was close to minimum ionising and considerably more uniform through the 
detector than positrons. 

Strict timing cuts were made to ensure that only events with single par- 
ticle interactions were selected, and that events with multiple overlapping 
events from the test beam were rejected. This was done since rate of events 
at the FD is so low (about 1 Hz of cosmic ray muons) there is little to be 
gained from comparing the response for multiple overlapping particle inter- 
actions occurring at the same time ("pile-up" events). Further timing cuts 
were made to ensure that the single particle interactions were measured by 
both systems, since the length of time the two systems were enabled was 
different. 

Positrons were selected by demanding the following: a TOF consistent 
with their velocity; a Cherenkov counter hit compatible with a single positron; 
coincidence in time of the TOF trigger and Cherenkov counter; and energy 
depositions only in the first 20 planes (to further reduce the number of pile- 
up events, specifically the case where an electron arrives in coincidence with 
a longer muon or pion). 

High energy muons were selected by demanding: a TOF consistent with 
their velocity; no signal in the Cherenkov counter; a particle track crossing all 
60 planes of the detector; and energy depositions consistent with a minimum 
ionising particle. 

4. Calibration 

The intrinsic calorimctric response of the ND and FD readout systems is 
different and thus it is necessary to calibrate their response independently. 
The readout system comparisons presented in this paper are dependent on the 
calibration and hence constitute a precise test of the calibration procedure's 
ability to remove the intrinsic response differences between the two readout 
systems. 
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The calibration of CalDet relies on an LED-based light injection sys- 



tem [13j and cosmic ray muons. A detailed description of the calibration 
techniques developed for MINOS using CalDet are given in [3]. A multi- 
stage procedure that converts the raw charge Qraw{,i,t,x) measured by strip 
i at time t and position along the length of the strip x into a corrected signal 
Qcor was used. Each calibration stage produced a numerical factor ("cali- 
bration constant"). Qcor is defined by the product of Qrawiht^x) and the 
calibration constant from each stage: 

Qcor Q raw (^^ ^ ^ ^ ^ -^(^5^) 

xU{i) X A{i,x) X S, 

where D, U, A and 5* refer to: 

Drift Correction D{i,t): To correct for temporal variations, data from 
the light injection system were used to track the response of each channel in 
the two readout system^. 

Uniformity Correction U{i): Cosmic ray muons that passed right 
through the detector were used to correct for differences in the overall light 
output of each scintillator strip. For the data presented herein, this calibra- 
tion was done once and a 2% statistical precision was achieved. 

Attenuation Correction A{i,x): A correction for light attenuation 
along the length of strips was made separately for the horizontal and vertical 
strips. This attenuation correction had to be made because the test beam 
particles hit the detector up to 2 cm away from the centre of the 1st plane 
depending on their momentum and particle type. Deviations from the centre 
of the detector meant that, without correction, the response measured by one 
readout system was systematically low or high with respect to the other. 

MC simulations of positrons striking the detector in a number of different 
positions were used to compute overall attenuation correction factors used 
for both data and MC. A correction factor of ±0.2%/cm for both horizontal 
and vertical strips was applied separately to the two readout systems to give 
the response as if the positrons struck the centre of the plane. 

Signal Scale Correction 5*: Stopping muons are used to set an overall 
energy scale [l^. A muon energy unit (MEU) is defined, which corresponds 
to the energy deposited by an approximately 0.8 GeV muon passing per- 



•^The light injection data were also used to determine the gain of each channel but this 
information was not directly used in determining Qcor- 
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pendicularly through a 1 cm thick MINOS scintillator plane (the average 
absolute energy deposition by such a muon is approximately 2 MeV). An 
energy deposition of 1 MEU yields about 4 PE per strip-end. The same scale 
factor is used for both readout systems; i.e. the MEU scale has no impact 
on the relative response comparison. However, the MEU scale is defined to 
allow calorimetric response comparisons among the three MINOS detectors. 

5. Calibration Uncertainties 

5.1. Drift Correction 

The drift in the absolute response of the two readout systems with time 
was of the order of several percent 0] over several weeks of running. This 
was mostly due to large temperature variation in the experimental area of up 
to 10°C. However, the absolute temperature dependence of the two readout 
systems was found to be almost identical and hence the error on the relative 
response was a negligible 0.1%. This was determined from 1 GeV/c test 
beam positron data sets taken at regular intervals. 

5.2. Uniformity Correction 

The performance of the uniformity correction is illustrated in Fig. [2] by 
plotting the average calibrated response per plane for test beam muons with 
momentum greater than 2 GeV/c. These muons travel through the detector 
by entering at the centre of the first plane and exiting within a 25 cm radius 
of the last plane's centre. The FD and ND readout systems have a 2.4% and 
3.0% spread in response from plane to plane respectively, which demonstrates 
the size of the systematic uncertainties on the uniformity correction. 

5.3. Non-Linearity Measurements 

The response of the readout system to light from the scintillator is not 
perfectly linear due to non-linearities in the PMT and electronics response. 
To incorporate the correct level of non-linearity into the MC simulation a data 
driven approach was taken. The light-injection system was used to measure 
the average non-linearity (1-measured/expected) of all readout channels for 
the two systems separately. These measurements were then used in the MC 
simulation to induce a non-linear response. 

The data from the light-injection system are shown in Fig. [3], where the 
non-linearity of each readout system is plotted as a function of the number 
of expected photoelectrons. PIN photodiodes that measured the quantity of 
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Figure 2: Average calibrated energy deposited versus plane by >2 GeV/c test beam muons. 
The error bars show the statistical error, thus the jitter from point to point is indicative of 
systematic error in the uniformity calibration. The observed response spread of the points 
is 2.4% and 3.0% for the FD and ND readout systems respectively. 
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Figure 3: Average non-linearity of the ND and FD readout systems versus the number of 
PE injected by the hght injection cahbration system. The difference in the non-Unearity 
between the two readout systems can be clearly seen. 
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injected light with a high degree of hnearity were used to determine the ex- 
pected number of photoelectrons. The ND PMTs are more non-hnear than 
the FD PMTs. A fit to the CalDet data gave -0.05%/PE and -0.03%/PE 
deviations from hnearity for the ND and FD PMTs respectively [sj. How- 
ever, as can be seen from Fig. [3], the PMT non-linearity does not explain all 
the data. The shape of the overall (PMT plus electronics) readout system 
non-linearity is quite different for the two cases. This occurs because of a 
saturation effect in the FD electronics that becomes significant for signals 
above about 100 PE. Beyond 200 PE the FD readout system can be seen to 
become significantly more non-linear than the ND system. Overall, the non- 
linearities of both readout systems are the same to within 2% up to about 
200 PE. 

No correction for non-linearity was applied to the data or the MC. This 
approach was taken in order to investigate and understand the relative effects 
of the non-linearity. Furthermore, by considering this worst case scenario it 
is possible to set a maximum error on the relative energy calibration between 
the two detectors. This maximum error can then be reduced by performing 
a non-linearity calibration in the Near and Far detectors. However, since the 
non-linearity of the two readout systems is similar (see above) , the systematic 
error on the relative energy scales is expected to be small for the range of 
energy depositions relevant to MINOS. 

6. Strip Occupancy and PMT Crosstalk Comparison 

The strips that register a response for a given particle interaction is not 
always the same for the two readout systems. Differences in strip occupancj0 
can arise due to PMT crosstalk, the average number of photons striking the 
PMTs and readout thresholds. The effect of crosstalk on the strip occupancy 
can be clearly demonstrated by using test beam muons since they only strike 
the central few strips all along the length of the detector. Fig. H] shows the 
relative strip occupancy for test beam muons with momentum greater than 
2 GeV/c. The central peak is mostly formed from strips that the muon 
actually passed through. In contrast, the peaks towards the edges of the 
detector (strip number 0-3 and 20-23) are dominated by crosstalk hits. The 
strip occupancy at the detector edges is different between the two readout 



^ Strip occupancy is a count of the number of times a specific strip is hit, normahsed 
by the total number of strips registering a hit. 
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systems due to the different cabling configurations and PMTs used. Overall, 
it can be seen that the strip occupancy is well reproduced by the MC for 
both readout systems. 

In addition to studying the strip occupancy, detailed studies were made 
of the crosstalk at the PMT level. The crosstalk arises primarily when pho- 
toelectrons injected on one pixel leak into the dynode chains of neighbour- 
ing pixels. Thus, crosstalk can be quantified as Q Neighbour I Q injected^ where 
Qinjected is the charge on the anode of the pixel where light was injected, 
and Q Neighbour is the charge that appears on the anodes of the neighbouring 
pixels. In-situ measurements of this crosstalk using muons in CalDet gave 
values of 5.2±0.1% and 4.8±0.1% for the FD and ND respectively 0. How- 
ever, since 16 of the 64 pixels on each ND PMT were not used, the observed 
crosstalk fraction associated with strips dropped to 3.2 ± 0.1%, consistent 
with expectation^. 

7. Calorimetric Response Comparison 

The calorimetric response of the readout systems are compared in two 
distinct ways: a "hit-by-hit" comparison and an "event-by-event" compari- 
son. In the hit-by-hit comparison the difference in the response from the two 
ends of individual strips is taken. For the event-by-event comparison, it is 
the sum total of the response of all the strip-ends read out by one system 
that is compared with the sum total of the response of all the strip-ends read 
out by the other system. In this paper, the differences are characterised by 
a relative response asymmetry, Aj\j-/jr. In general form, Aj^/jr is written as 



where Af and JF are the calorimetric response of the ND and FD readout 
systems respectively. The data are presented here as the average hit-by- 
hit asymmetry, < A^^jj^ >, and as the average event-by-event asymmetry, 
< A^J'^'^ >. It should be noted that in order for a comparison to be made 
in the hit-by-hit case, a response above threshold on both ends of the strip 



^The naive expectation is for the crosstalk associated with strips to drop by 16/64 to 
3.6%. However, effects arise from exactly which 48 of the 64 pixels are connected to strips 
since some pixels receive more crosstalk than others. 
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Figure 4: Transverse hit profile of >2 GeV/c test beam muons as measured witli tlie ND 
(top) and FD (bottom) readout systems for data and MC. The central peak is mostly 
formed from strips that the muon actually passed through. Whereas, in contrast, the 
peaks towards the edges of the detector (strip number 0-3 and 20-23) are dominated by 
crosstalk hits. These data are from planes with strips that run horizontally. 
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was required. In contrast, strips with a response at just one end were in- 
cluded in the event-by-event comparison. Depending on the particle energy, 
approximately 90-95% of the total event energy comes from strips with a 
hit on both ends. The remainder comes from strips with a response on only 
one end and that fraction is consistent with PMT crosstalk, fluctuations in 
the production of photoelectrons, and readout thresholds. 

7.1. Hit-by-hit Response Comparison 

The average hit-by-hit asymmetry, <Aj^Jyr>, is shown in Fig. [5] as a func- 
tion of energy deposition. The data from positrons of between 0.6-6.0 GeV/ c 
were used to cover the range of energy depositions from 0-36 MEU. The 
uniformity correction forces the asymmetry to be zero at around 1.5 MEU, 
which is the average energy deposition per strip of the cosmic ray muons 
used for the calibration of CalDet. Variation of < > as a function of 

energy deposition is best described by considering 2 distinct regions: above 
and below 4 MEU. 

<A^Jyr> changes by up to 2% below 4 MEU with data and MC agreeing 
to 1%. The general trend is that < A^^jj^ > increases by about 1.5% be- 
tween 0.5-4 MEUs of energy deposited, which means that the ND response 
is increasing faster with energy deposited than the FD response. This vari- 
ation is caused by the difference in the way the readout threshold is applied 
for the two readout systems: the ND electronics applies a threshold to each 
18.8 ns integration sample separately, whereas the FD electronics only applies 
a threshold to the total integration sample (see Section [2] for more details). 
As the energy deposited increases, the fraction below readout threshold tends 
towards zero and so this effect is no longer significant. This trend is well de- 
scribed by the MC simulation, which shows a maximal deviation of less than 
0.6% from the data in all but the lowest energy bin. 

Between 4-15 MEU the slope goes the opposite way with < A^^j^^ > 
decreasing with energy deposited by 2.5%. This is caused by the fact that 
the ND PMTs are less linear than the FD PMTs, as shown in Fig. [31 For 
energy depositions between 18-36 MEU, the trend is reversed with the slope 
of < Aj^Jjr > becoming positive again. The location of this inflection is 
consistent with the light level at which the FD readout system non-linearity 
starts to increase faster than that of the ND (see Fig. [3]). The data and 
MC agreement in the region between 4-36 MEU is better than 1%. It is 
important to note that the results in this region quantify the maximum error 
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Figure 5: Average hit-by-hit asymmetry ratio, < ^_^/jr >, as a function of the average 
energy deposition in the scintiUator strip hit for data and MC. It can be seen that difference 
between the response of ND and FD readout systems changes by up to 2.5% depending 
on the energy deposition. 



16 



in the case when a non-hnearity correction was not apphed. Whereas, when 
a non-hnearity correction is apphed the differences between the two readout 
systems at energy depositions above 4 MEU will be reduced. 

In summary, the hit-by-hit data presented here show energy dependent 
differences between the calorimetric response of the two readout systems. 
These differences are shown to be caused by threshold effects as well as 
PMT and electronics non-linearity. The fact that the differences in the two 
readout systems are well simulated is particularly important for the relative 
calibration between the MINOS detectors. This result shows that as long 
as data and MC are made to agree at the energy of the muons used for 
calibration (around 1.5 MEU) then the MC can be used to correct for the 
relative differences between the ND and FD readout systems. 

1.2. Event-by-Event Response Comparison 

To make the comparison of the event-by-event response of the two readout 
systems, data from test beam positrons between 0.6-6.0 GeV/c were used. 
Fig. [6] shows the average event-by-event asymmetry, <A^I^^^>, as a function 
of the test beam momentum. It can be seen that there is no significant 
energy dependence. This is an important demonstration that the differences 
between the two readout systems seen in the hit-by-hit comparison do not 
introduce an overall relative energy dependent bias between the ND and FD 
(despite the absence of a non-linearity correction in the data presented here). 
It should be noted that a correction of 2.0% for the difference in crosstalk 
between the two readout systems was made in data and MC (see Section [6]). 

The average value of <y4_^"^*> between 0.6-6.0 GeV/c is 

1.2 ± 1.3% for data, and 
0.3 ± 0.3% for MC. 

The statistical errors on both data and MC are negligible. The systematic 
uncertainties arise from the following sources: the uncertainty in the test 
beam position used to make the attenuation correction; the correction for 
differing quantities of crosstalk associated with strips; and bias in the uni- 
formity correction. 

The uncertainty on the position of the centre of the test beam in both 
data and MC gives a 0.25% error on < A^J^^'^ >. The uncertainty on the 
difference in the cross talk between the two readout systems gives a 0.14% 
error. For MC the systematic error from the uniformity correction is by 
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Figure 6: Average event-by-event asymmetry ratio, < >, as a function of the test 

beam momentum. The error bars show the total systematic error (statistical errors are 
negligible) . 
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definition zero since the same calibration constants are used to both induce 
differences between the response of the strips in the simulation and then to 
calibrate out those differences. 

For data the systematic uncertainty is dominated by the uniformity cal- 
ibration. In order to estimate the systematic uncertainty on the unifor- 
mity calibration muons from the test beam with a momentum of more than 
2 GeV/c were used. The value of < A^J^^ > obtained using these muons is 
shown as a function of position along the detector in Fig. [71 The data show 
a 1.1% difference between the first third of the detector and the back two 
thirds. In contrast, the asymmetry in the MC simulation is close to uniform 
throughout as expected. The presence of the asymmetry at the front of the 
detector in data is attributed to a bias induced by the uniformity correc- 
tion. The measured asymmetry value of 1.1% is used as an estimate of the 
systematic error on the uniformity calibration and is the dominant error on 
<A^J^>. 

7.3. Readout System Resolution 

The data described in this paper were obtained with the scintillator strips 
instrumented at one end with the ND readout system and at the other end by 
the FD system. The advantage of this approach was that it made the com- 
parison of the two systems insensitive to fluctuations in the energy loss of the 
particles themselves and the simulation of such fluctuations in the MC. The 
fluctuations observed in A^JJ^^ from event to event are due to fluctuations in 
the detector readout and the absorption/re-emission of photons in the WLS 
fibre. The dominant source of fluctuations arises in the conversion of photons 
to photoelectrons at the PMT face. For example, a 1 GeV/c positron yields 
about 180 PE and leads to a 15% width. The width of the distribution of 
A^JJ^^ is shown in Fig. [8] as a function of positron momentum. It can be seen 
that the size of the fluctuations decreases with energy and is well reproduced 
by the MC simulation. 

8. Conclusions 

The MINOS Calibration detector has acquired test beam data in a config- 
uration that allowed direct comparison of the two different readout systems 
used for the Near and Far detectors. These data were used to understand the 
systematic differences between the readout systems and to demonstrate that 
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Figure 7: < ^_^/*jr > computed with test beam through going muons along the detector 
axis (statistical errors only). Data exhibits a difference in < A^Jy^ > between the first 
third and the rest of the detector, which is consistent with calibration systematic bias. 
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Figure 8: The readout system resolution (width of A^^^p*) as a function of positron 

momentum. This figure demonstrates the size of the fluctuations that occur in the readout 
when the same physics event is observed by both Near and Far systems. It can be seen 
that the size of the fluctuations decreases with energy and is well reproduced by the MC 
simulation. 
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the differences can be controlled to the level required by MINOS for neutrino 
oscillation measurements. 

The results for the hit-by-hit calorimetric response comparison can be 
divided into two regions: energy depositions above and below 4 MEU. The 
relative hit-by-hit response of the two readout systems below 4 MEU changes 
by up to 2% with the MC reproducing the trend to 1%. The cause for 
this difference is demonstrated to be dominated by threshold effects in the 
electronics. This affect could potentially bias the relative energy calibration 
in MINOS. One way to avoid this bias is to calibrate the MC so that it agrees 
with the data in each detector separately. Above about 4 MEU the hit-by-hit 
response of the two systems differs by up to 2.5% with the MC reproducing 
the trends to better than 1%. This 2.5% variation is caused by differences in 
the non-linearity of the two types of PMT and electronics used for the Near 
and Far detectors. 

The event-by-event comparison showed that the differences in the overall 
calorimetric response of the two readout systems were consistent with zero to 
an accuracy of 1.3% and 0.3% for data and MC respectively. Furthermore, 
no significant energy dependence was observed in positron data sets taken 
with test beam momentum settings between 0.6-6 GeV/c. 

Differences in strip occupancy were shown to arise due to effects of crosstalk 
between pixels on the two types of multi-anode photomultiplier tubes used. 
This crosstalk, which is different for the two readout systems was demon- 
strated to be well reproduced by the MC simulation: both in the fraction 
of the total energy deposition dispersed as crosstalk and in its distribution 
across the pixels. 

In the MINOS Near and Far detectors it is neutrino interactions that are 
reconstructed and that involves a multiplicity of particles in the hadronic 
shower. However, it is the combination of the result presented here for single 
positrons and the demonstrated accuracy of the MC simulation that is most 
important. These two results combined give confidence that any effects aris- 
ing from the two different readout systems will be small enough so that they 
do not significantly impact the neutrino oscillation measurements made by 
MINOS. 
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